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Erneux,  Mixed  state  effects  in  waveguide  electro-absorbers  based  on  quantum 
dots,  Appl.  Phys.  Letters  99,  171103  (2011) 

2.  Christian  Otto,  Kathy  Ltidge,  Evgeniy  Viktorov  and  Thomas  Erneux, 
Quantum  dot  laser  tolerance  to  optical  feedback,  In  “Nonlinear  Laser  Dy¬ 
namics”,  K.  Ltidge  (Ed),  Wiley-VCH  Weinheim,  Germany  (2012) 


In  (1),  multi-pulse  heterodyne  pump-probe  measurements  are  used  to  in¬ 
vestigate  the  reverse  bias  dynamics  of  InAs/GaAs  quantum  dots  in  a  waveguide 
structure.  It  is  the  last  publication  of  a  series  of  pump-probe  experiments 
done  in  collaboration  with  the  group  of  G.  Huyet  at  Cork  [1]- [4] . 

Publication  (2)  is  a  chapter  of  a  book  collecting  different  contributions 
on  Laser  Dynamics.  It  considers  the  stability  condition  derived  by  Mork  et 
al.  [5]  given  by 


k  <  kc  = 


v7!  +  a2 


(1) 


where  a  is  the  linewidth  enhancement  factor  and  T®w  is  defined  as  the  damp¬ 
ing  rate  of  the  relaxation  oscillations  multiplied  by  the  diode  cavity  round- 
trip  time.  Eq.  (1)  was  previously  suggested  by  Helms  and  Petermann  [6] 
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as  a  simple  analytical  criteria  for  tolerance  with  respect  to  optical  feedback. 
Helms  and  Petermann  [6]  also  proposed  the  empirical  law  given  by 


=  pQW  +  q2 

n/2 


(2) 


Both  Eqs.  (1)  and  (2)  are  used  in  current  experimental  studies  of  quantum 
dot  (QD)  lasers  subject  to  optical  feedback.  This  has  motivated  the  analysis 
of  two  different  rate  equation  models  for  QD  lasers.  In  both  cases,  we  derive 
an  expression  that  has  the  same  format  as  1  but  with  a  expression  for  the 
damping  rate  YqD  >  TqW. 


A  combined  experimental-theoretical  paper  describing  the  response  of  QD 
lasers  in  turn-on  experiments  has  been  submitted.  Another  paper  on  optical 
feedback  is  in  preparation. 
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A.  Viktorov,  T.  Erneux,  and  S.  P.  Hegarty,  Phys.  Rev.  E84,  026208  (2011) 

Paper  (1)  reports  on  an  experimental  and  theoretical  study  of  a  single 
mode  quantum  dot  (QD)  laser  subject  to  injection.  The  stability  diagram 
considerably  differs  from  the  one  known  for  the  quantum  well  (QW)  semi¬ 
conductor.  Experimentally,  two  features  stand  out;  the  first  is  an  absence  of 
instabilities  resulting  from  relaxation  oscillations,  and  the  second  is  the  ob¬ 
servation  of  a  region  of  bistability  between  two  locked  solutions.  Using  rate 
equations  appropriate  for  quantum-dot  lasers,  we  analytically  determine  the 
stability  diagram  in  terms  of  the  injection  rate  and  frequency  detuning. 

At  the  TU-Berlin,  a  microscopically  based  rate  equation  model  for  QD 
lasers  is  studied  mainly  numerically.  The  model  separately  treats  the  dy¬ 
namics  of  electrons  and  holes,  and  the  carrier-carrier  scattering  rates  depend 
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nonlinearly  on  the  wetting  layer  carrier  densities.  Paper  (2)  summarizes  more 
than  two  year  efforts  to  simplify  the  five  rate  equations  in  a  form  for  which 
analytical  expressions  can  be  obtained  for  the  relaxation  oscillation  frequency 
and  damping  rate.  . 

Paper  (3)  describes  new  experimental  observations  of  noise-induced  ex¬ 
citable  pulses  for  an  injected  QD  laser.  The  authors  use  a  third  order  phase 
equation  instead  of  Adler  equation  for  the  interpretation  of  their  results. 
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Experiments  on  a  quantum  dot  (QD)  DFB  laser  subject  to  optical  injec¬ 
tion  have  been  realized.  The  stability  diagram  has  been  determined  in  terms 
of  the  slave-master  frequency  detuning  and  the  optical  injection  rate.  The 
domains  of  stability  strongly  differ  from  the  same  domains  observed  for  the 
conventional  quantum  well  laser.  We  have  considered  a  three-variable  model 
for  the  QD  laser  [1]  in  the  case  where  it  admits  a  strong  damping  of  the 
relaxation  oscillations.  We  explain  numerically  and  analytically  the  unusual 
stability  properties  observed  experimentally.  A  publication  with  the  group 
of  G.  Huyet  (Cork)  in  Optics  Letters  summarizes  the  results  and  acknowl¬ 
edge  the  support  of  the  AF  [2] .  Another  manuscript  has  been  submitted  for 
publication. 

Parallel  to  the  collaboration  with  Cork,  work  has  been  initiated  with  the 
group  of  E.  Scholl  (TU  Berlin)  in  order  to  simplify  a  five  variable  model  used 
by  the  group  [3].  Our  objective  is  to  determine  if  there  exist  similar  features 
between  this  model  and  the  three  variable  model  used  in  Cork.  A  manuscript 
summarizing  the  asymptotic  analysis  is  in  preparation. 

By  the  end  of  2009,  the  group  in  Brussels  became  partner  of  the  French 
ANR  project  called  TELDOT  devoted  to  the  development  of  quantum  dot 
lasers  for  Telecom  applications.  To  this  end,  work  has  been  started  with 
the  group  of  P.  Besnard  (Lannion).  This  group  is  now  performing  injection 
experiments  using  a  quantum  dash  DFB  laser  fabricated  by  Alcatel.  This 
study  could  be  particularly  relevant  for  the  AF  which  has  recently  supported 
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experiments  using  a  quantum  dash  laser  in  a  Fabry- Perot  cavity  [4] . 
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Abstract 

Recent  experiments  have  explored  the  response  of  a  single  mode 
optically  injected  QD  laser.  The  experiments  identify  a  domain  of 
bistability  between  two  steady  states  at  negative  detuning  and  a  Hopf 
bifurcation  close  to  the  locking  threshold  at  positive  detuning.  More¬ 
over,  the  laser  is  stable  for  all  injection  levels  at  zero  detuning.  These 
observations  lead  to  a  stability  diagram  injection  amplitude  versus  de¬ 
tuning  that  contrasts  to  the  diagram  of  a  quantum  well  semiconductor 
laser.  The  bifurcations  responsible  for  the  QD  laser  stability  diagram 
are  explained  analytically  by  using  a  model  appropriate  for  QD  lasers. 


1  Introduction 

Semiconductor  lasers  (SLs)  have  become  the  optical  source  of  choice  in  many 
applications  due  to  their  high  efficiency,  simplicity  of  modulation,  and  small 
size.  However,  in  some  applications  where  low-intensity  noise  level  are  re¬ 
quired,  they  suffer  from  wide-bandwidth  intensity  fluctuations  that  are  en¬ 
hanced  by  their  inherent  relaxation  oscillations  (ROs).  The  phenomenon  of 
ROs  is  familiar  in  the  laser  physics  community.  When  a  laser  is  perturbed 
from  its  steady  state  operation,  it  does  not  immediately  return  to  its  original 
position.  Either  the  laser  quickly  approaches  equilibrium  or  it  slowly  decays 
to  its  stable  steady  state  like  a  damped  oscillator.  These  lasers  are  labelled  as 
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Class  A  and  Class  B,  respectively  [1].  Class  A  lasers  include  the  Ar,  He-Ne, 
and  dye  lasers  while  Class  B  lasers  include  most  of  the  lasers  used  today  such 
as  CO2,  solid-state,  and  semiconductor  lasers.  Subject  to  optical  injection 
Class  A  and  Class  B  lasers  exhibit  quite  different  stability  properties.  Class 
B  lasers  admit  a  rich  number  of  sustained  pulsating  intensity  regimes  which 
have  been  studied  systematically  during  the  last  decade  for  semiconductor 
and  solid  state  lasers  (see  [2]  for  a  recent  review).  Class  A  lasers  free  of  ROs 
are  much  more  stable  [3].  Recent  efforts  have  concentrated  on  increasing 
the  photon  lifetime  above  the  carrier  lifetime  to  suppress  ROs  in  SLs.  This 
can  be  achieved  by  either  increasing  the  cavity  length  or  the  cavity  finesse. 
The  first  technique  has  been  successfully  applied  with  a  several-meter-long 
cavity  in  semiconductor  laser  [4] .  The  second  alternative  is  technically  easier 
and  has  been  explored  using  half  of  a  vertical  cavity  surface  emitting  laser 
(|— VCSEL)  in  a  short  external  cavity  [5,  6]. 

We  consider  optical  injection  of  a  single  mode  distributed  feedback  (DFB) 
quantum  dot  laser  (QDL).  These  lasers  have  become  increasingly  pervasive 
in  recent  years  and  studies  have  already  revealed  some  dynamical  properties 
which  demonstrate  their  superiority  for  applications  [7] .  A  particular  feature 
of  these  devices  is  an  unusually  high  damping  of  the  relaxation  oscillations 
(ROs)  [8]  in  comparison  to  their  bulk  and  quantum  well  (QW)  counterparts. 
Analytical  studies  have  shown  that  the  damping  rate  and  the  RO  frequency 
can  have  the  same  timescale  and  may  even  be  overdamped  [9].  This  high 
damping  has  been  cited  as  the  principal  reason  for  the  increased  stability 
of  such  devices  subject  to  optical  feedback  [10],  optical  injection  [11],  and 
mutual  coupling  [12]  configurations.  We  determine  an  experimental  stability 
diagram  and  note  that  it  is  considerably  different  to  that  of  a  conventional 
quantum  well  laser  showing  strong  similarities  to  the  stability  diagram  of 
a  Class  A  laser.  In  particular,  the  injected  QDL  exhibits  stability  for  ar¬ 
bitrary  value  of  the  injection  rate  provided  the  detuning  is  sufficiently  low 
and  bistabiliy  between  two  coexisting  steady  states  at  negative  detuning. 
These  observations  are  substantiated  analytically  by  studying  rate  equations 
appropriate  for  a  QD  laser. 


2  Experiments 

The  slave  DFB  QDL  used  was  a  five  layer  structure  grown  by  solid  source 
MBE.  It  consisted  of  2.4  InAs  monolayers  topped  with  5  nm  GalnAs,  stacked 
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in  a  400  nm  thick  optical  cavity.  A  35nm  GaAs  spacer  is  used  between 
the  QD  layers.  Optical  confinement  was  achieved  using  AlGaAs  cladding 
layers.  The  experimental  setup  is  similar  to  the  one  described  in  [13].  Of 
particular  interest  is  the  large  area  of  stable  locked  operation.  Regions  of 
bistability  between  locked  and  periodic  operation  exist  for  both  negative  and 
positive  detunings  while  for  negative  detuning,  there  is  a  domain  of  bistable 
operation  in  which  two  locked  steady  states  coexist  and  the  laser  displays 
a  noise  induced  switching.  The  transition  between  the  two  steady  states  is 
sharp  and  the  relaxation  dynamics  includes  only  one  spiking  oscillation.  It 
suggests  that  the  decay  of  the  relaxation  dynamics  occurs  at  the  same  time 
scale  as  that  of  the  RO  frequency  which  is  typically  the  case  for  a  Class  A 
laser. 

There  are  a  number  of  fundamental  differences  between  our  stability  di¬ 
agram  and  the  diagram  reported  in  [14]  for  the  conventional  DFB  QW  laser 
(reviewed  in  detail  in  [2]).  In  the  injected  QW  laser,  there  is  only  one  Hopf 
bifurcation  line  that  crosses  twice  the  zero  detuning  line.  Here  we  found  two 
Hopf  bifucation  lines  that  never  cross  the  sero  detuning  line.  Moreover,  the 
coexistence  of  two  stable  locking  states  for  our  QW  laser  is  not  possible  for 
a  QW  laser  except  close  to  the  laser  threshold.  These  differences  suggest  a 
significant  impact  from  the  nonlinear  capture  dynamics  in  QDLs  provided  by 
the  Pauli  blocking  factor  which  do  not  exist  in  the  conventional  formulation 
of  a  class  B  QW  laser.  Instead,  we  note  a  similarity  between  our  stability 
diagram  and  the  one  for  a  class  A  laser  [3].  It  suggests  that  QDLs  may 
exhibit  both  class  A  and  class  B  dynamics  depending  on  the  carrier  capture 
parameters  as  shown  in  [9]  by  exploring  different  asymptotic  limits  of  a  three 
variable  rate  equation  model. 


3  Theory 


Our  rate  equations  for  a  QD  laser  subject  to  an  injected  signal  consist  of 
three  equations  for  the  complex  electric  field  E,  the  occupation  probability 
in  a  dot  p,  and  the  carrier  density  n  in  the  wetting  layers,  scaled  to  the  QD 
density.  They  are  given  by  [11] 


E' 

P 


-(1  +  ioP)  [— 1  +  g{2p  —  1)]  E  +  T exp(iAf), 


=  V 


Bp(l-p)-p-(2p-l)\Ef 


(1) 

(2) 
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2Bn(l  —  p)\ . 


(3) 


n  =  77  [  J  —  n  — 

Prime  means  differentiation  with  respect  to  t  =  t' /rph  where  rph  is  the  photon 
lifetime.  The  factor  1  —  p  is  the  Pauli  blocking  factor.  The  factor  2  in  Eq. 
(3)  accounts  for  the  spin  degeneracy  in  the  quantum  dot  energy  levels.  J  is 
the  pump  current  per  dot  and  a  is  the  linewidth  enhancement  factor.  The 
control  parameters  are  the  frequency  detuning  A  defined  as  the  frequency  of 
the  master  laser  minus  that  of  the  slave  laser  and  the  injection  rate  F.  The 
fixed  parameters  B  and  77  are  ratio  of  basic  time  scales  and  are  defined  as 
B  =  TT~aP  and  77  =  t^t-1  where  r  and  Tcap  denote  the  carrier  recombination 
and  capture  times,  respectively.  Typical  values  are  r  =  1  ns,  rph  =  2  ps,  and 
rcap  =  10  ps  which  imply  B  =  102  and  77  =  2  x  10-3. 

As  suggested  by  Goulding  et  al  [11],  we  shall  consider  the  value  g  =  1.01 
for  which  a  good  agreement  between  theory  and  experiments  is  observed. 
In  the  case  of  the  solitary  laser,  the  product  B(g  —  1)  appears  in  both  the 
steady  state  expressions  and  in  the  characteristic  equation  [9].  Therefore,  we 
need  to  take  into  account  the  relative  values  of  B  and  g  —  1.  Specifically,  we 
propose  an  asymptotic  analysis  of  Eqs.  (l)-(3)  valid  in  the  limit  e  =  g  —  1 
small  keeping  Be  as  an  0(1)  quantity.  After  introducing  g  =  l  +  e  into  Eq. 
(1),  the  expression  in  brackets  becomes  [—2  +  2p  +  e(2 p  —  1)]  and  suggests 
to  introduce  p  =  1  +  eu  in  order  to  balance  all  terms.  The  expression  in 
brackets  then  is  proportional  to  e  which  motivates  introducing  the  slow  time 
scale  s  =  et.  From  Eqs.  (l)-(3),  we  obtain  the  following  equations  for  E.  u. 
and  n 


E' 

=  ^(1  +  ia)  [1  +  2m(1  +  e)]  E  +  7 exp(i<5s), 

(4) 

v! 

=  e~2g  —Benu—l  —  eu  —  (l  +  2ue)\E\2  , 

(5) 

n' 

=  £-1?7  [  J  —  n  —  2Benu\ 

(6) 

where  prime  now  means  differentiation  with  respect  to  s.  The  control  pa¬ 
rameters  are  7  =  £_1T  and  8  =  e_1A.  Since  e~2  £-1  as  e  — >  0,  u  is 

faster  than  n  and  we  eliminate  u  adiabatically  from  Eq.(5).  Specifically, 
we  find  u  =  —(1  +  E2) /(Ben)  as  e  — ■>  0.  Introducing  the  decomposition 
E  =  R exp(ibs  +  id),  Eqs.  (4)-(6)  reduce  to  the  following  equations  for  R,  0, 
and  n 

R'  =  \  1~2ilBef)  B  +  TcosW-  (7> 
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(9) 


If  7  =  0,  we  find  that  the  laser  threshold  of  the  solitary  laser  appears  at 
J  =  Jth  where 

Jtfc  =  2  +  2/(£e).  (10) 

We  next  assume  that  J  >  Jth-  From  the  steady  state  equations,  we  determine 
R2  =  R2(y)  as  (in  implicit  form) 


7 2  =  R2  F2  +  (-8  +  aFf 


(11) 


where 

F  =  l(! _ 2(!  +  R2) _ 

2\  Be(J-  2(1 +  fl2)) 

From  the  linearized  equations,  we  then  obtain  the  characteristic  equation 
for  the  growth  rate  A.  It  is  given  by 


A3  -|-  cqA  T  0,2^  F  cis  —  0, 

where 

(i\  =  — G  —  F  F  e  ^  7 , 
a2  =  GF  F  (F  -  G)a{8  -  aF)  F  (8  -  aF)2 
-£~1r/(G  F  F)  F  e-'viF  -  G)(l  -  2 F), 
=  GF  F  (F  -  G)a{6  -  aF)  F  (8  -  aF)2 
F(F  -G)(  1  -  2 F)  (a{6  -  aF)  -  F) 

and 

If  2J1F3R2)  \ 

2\  Be  ( J  —  2(1  +  R?)  J 


(13) 

(14) 

(15) 

(16) 

(17) 


The  Routh-Hurtwitz  stability  conditions  are  ai«2  —  c/3  >  0,  «i  >  0,  and 
«3  >  0.  The  saddle-node  bifurcation  point  satisfies  the  condition  03  =  0.  The 
Hopf  bifurcation  point  satisfies  the  condition  Gqa2  —  03  =  0.  Both  conditions 
are  quadratic  equations  in  (8  —  aF).  We  first  determine  8  —  aF  as  a  function 
of  R2  and  then  8  as  a  function  of  7,  using  (11).  The  stability  diagram  is 
shown  in  Fig.  1.  Only  the  Hopf  bifurcation  points  from  a  stable  steady  state 
are  shown  (0702  —  03  =  0,  oq  >  0,  and  03  >  0). 
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Figure  1:  Stability  diagram.  SN  and  H  denote  the  saddle-node  and  Hopf 
bifurcation  points,  respectively.  The  shaded  region  corresponds  to  the  do¬ 
main  of  steady  state  bistability.  The  values  of  the  parameters  are  g  =  1.01, 
B  =  102,  r/  =  2  x  10-3,  a  =  1.2,  and  J  =  1.2  Jth  =  4.8  ( Jth  =  4).  The  dots  are 
Fold- Hopf  points  where  Hopf  and  SN  bifurcation  lines  merge.  Inset:  stability 
diagram  for  a  injected  Class  A  laser  (Eq.  (1)  in  [3]  with  Tk-1  =  1.2,  a  =  1.2, 
(3  =  1,  T  — ►  ctSk-1,  A  — ■>  —  Afhc-1,  and  t  — >  nt). 


6 


Distribution  A:  Approved  for  public  release;  distribution  is  unlimited. 


A 


A 


Figure  2:  Bifurcation  diagrams  of  the  stable  steady  and  time-periodic  solu¬ 
tions.  The  extrema  of  R  are  shown  as  functions  of  the  detuning  A.  The 
complete  S-shaped  branch  of  steady  states  is  shown  by  a  broken  line.  Same 
values  of  the  parameters  as  in  the  previous  figure,  (a)  partial  bistability  for 
T  =  0.0012  and  (b)  full  bistability  for  T  =  0.002.  The  figures  have  ben 
obtained  by  scanning  the  detuning  back  and  forth. 
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The  stability  diagram  in  Fig.  1  is  qualitatively  similar  to  the  experi¬ 
mental  mapping.  Both  the  experimental  and  analytical  stability  diagrams 
predict  stable  locking  for  arbitrary  values  of  the  injection  rate  provided  the 
detuning  is  sufficiently  small.  The  domain  of  multiple  steady  state  dominates 
for  negative  detunings  (bounded  by  the  lines  SNi  and  SNi  in  Figure  1).  At 
positive  detuning,  a  single  Hopf  bifurcation  (H2  in  Fig.  1)  emerges  from  a 
single  steady  state.  The  absence  of  Hopf  bifurcations  at  low  injection  levels 
has  previously  been  noted  allowing  the  observation  of  excitable  pulses  for 
both  negative  and  positive  detunings  [13] .  [A  domain  of  bistability  between 
locked  states  and  periodic  solutions  for  lower  injections  is  predicted].  At 
negative  detuning  and  for  relatively  high  injection  rates,  there  is  a  domain  of 
bistability  between  two  locked  states.  This  phenomenon  is  not  possible  for 
QW  SLs  except  if  the  pump  current  is  very  close  to  its  threshold  value[15]. 
The  bistability  phenomenon  is  here  possible  because  of  a  Hopf  bifurcation 
that  stabilizes  the  lower  intensity  branch  (Hi  in  Fig.  1).  Two  bifurcation 
diagrams  of  the  stable  steady  and  periodic  regimes  are  shown  in  Figure  2  il¬ 
lustrating  the  case  of  partial  and  full  bistability.  They  have  been  determined 
numerically  from  Eqs.  (7)-(9). 

The  Hopf  bifurcation  curves  do  not  cross  the  zero  detuning  line  as  for  QW 
SLs.  Consequently,  our  QD  SL  exhibits  higher  stability  properties  which  are 
important  for  some  practical  applications.  We  should  however  emphasize 
that  our  QD  laser  is  equivalent  to  a  Class  A  laser  because  it  verifies  the 
scaling  law  (g  —  1)B  =  0(1).  Other  scalings  are  possible  because  of  the  large 
diversity  of  QD  structures  that  are  currently  designed  possibly  leading  to 
different  conclusions  concerning  the  RO  damping  and  the  effects  of  injections. 
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